Robustness of regularities for energy centroids in the presence of random interactions 
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In this paper we study energy centroids such as those with fixed spin and isospin, those with 
fixed irreducible representations for both bosons and fermions, in the presence of random two-body 
and/or three-body interactions. Our results show that regularities of energy centroids of fixed spin 
states reported in earlier works are very robust in these more complicated cases. We suggest that 
these behaviors might be intrinsic features of quantum many-body systems interacting by random 
forces. 
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In 1998 Johnson, Bertsch and Dean obtained a prepon- 
derance of spinP^'''*'' = 0+ ground states for even-even 
nuclei in the presence of random two-body interactions 
0. Since then, there have been a lot of efforts towards 
understanding this observation. Studies along this line 
were reviewed in Refs. I^l^. 

Although recently there were progresses 0; IE in 
evaluation of ground state energies, finding the ground 
state by a simple approach is usually difficult in the pres- 
ence of random interactions. Thus one may study other 
quantities which are relatively simple. Along this line, 
much attention [7-13] has been paid to spin / energy 
centroids (defined by the average energy of spin I states 
and denoted by Ej). Main results are reviewed briefly as 
below. 

(1) From numerical experiments by using the TBRE, 
it was found in Ref. |1] that i?/'s with / ~ I^un or 
/ ~ /max have large probabilities to be the lowest while 
those with other / have very small probabilities to be the 
lowest. Roughly speaking, there are ~ 50% of the cases 
for which Ei with / ~ /min [I — /max) is the lowest. We 
define (i?/)min ((£'/)niax) as the value obtained by aver- 
aging El over the subset where Ei^i^.^ {Ei^i^,^^) is the 
lowest energy. Ref. 8j also demonstrated that (i?/)inin — 
CI{I + 1) and {E}) 

max — c\ (/max + 1) -/(/+!)], 

where the value of coefficient C depends on the active 
single-particle orbits and the choice of the ensemble. We 
have studied cases of single-j shell configurations as well 
as many-j shell configurations in which shells are de- 



noted by ji, j2, • • •, etc. For the TBRE, C ~ l/i'iY.ijf)- 
The regularities of Ej were argued in Ref. Q by as- 
suming that two-body coefficients of fractional parent- 
age (cfp's) behave like randomly, and the behavior of 
{Ei)niin — C'I{I -\- 1) was reproduced for four fermions in 
a. j — 17/2 shell under this simple assumption. 

The above regularities of Ei which are stable for single- 
closed shell (both single-j and many-j shells) were found 
in Ref. Q to be robust even for many-j shells where each 
orbit can have positive or negative parity and for systems 
with isospins. For cases of nucleons in many-j shells (dif- 
ferent j can have different parity), the value of coefficient 
in the relation (i?/±)min — C^/'''(/^-|-l) are sensitive 
to the j values but not to parity or isospin. Very recently, 
one of the authors of the present paper, Kota, studied in 
Ref. |lC| energy centroids with fixed irreducible repre- 
sentations of some of the group symmetries of the inter- 
acting boson models |0, [3| such as the sd interacting 
boson model (IBM), the sd IBM with isospin, etc. It was 
found that the lowest and highest irreducible representa- 
tions carry most of the probability for the corresponding 
centroids to be the lowest in energy. This is a general- 
ization of results found numerically in Refs. Q • 

(2) Mulhall, Volya, and Zelevinsky assumed in Ref. 

the geometric chaos (quasi-randomness in the process 
of angular momentum couplings) and derived a linear 
relation between Ej and /(/ -I- 1); The same result was 
derived in Ref. by resorting to the group structure 
of J7(2j + 1) D 0(3) for n fermions in a single-j shell. Let 
us define single-j Hamiltonian 



H = 



V2J+T 



(0) 
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The formula of Ej of Refs. [lllIT^ was written as follows, 
:e7 = ^(2J+1)G'j 



2j + l 



2j-20- + l)2(2j + l)2 



0(/2(/ + l) 



(2) 



where 0{P{I +1)'^) refers to higher / terms which seem 
negligible. The first term of this formula is a con- 
stant independent of /. The second term of Eq. (O 
is proportional to /(/ + 1). Thus we have the relation 
Ej ^ Eq + CI{I +1). However, the value of coefficient C 
thus obtained was found in Ref. to be systematically 
smaller than those obtained by the TBRE or empirical 
formula C ~ l/(4j^); and furthermore, even for systems 
in which one cannot assume randomness of the geomet- 
ric chaos or randomness of the cfp's, a similar pattern 
was found to occur. Therefore, the arguments of Refs. 
[E^ill^ are just part of the full story and a sound un- 
derstanding is not yet available. It is interesting to note 
that one can obtain the /(/ -I- 1) term of Eq. (j^J with an 
additional factor of ^ when one transforms the single-j 
hamiltonian of Eq. into its particle-hole form via the 
Pandya transformation. 

The purpose of this paper is as follows. First, we dis- 
cuss energy centroids with fixed spin / and isospin T, 
denoted by Ej^ts. Although Ref. Q discussed S/'s for 
systems with isospin, i?/'s with different isospin T were 
mixed. Second, we study a Hamiltonian with random 
three-body interactions for sd bosons, while earlier works 
studied the property of i?/'s by using random two-body 
interactions. Third, we study with a fixed irre- 

ducible representation {/}, as an extension of the work 
in Ref. [T^ . These results are discussed by using propa- 
gation equations. 

This paper is organized as follows. In Sec. II, we 
discuss results of iJ/.r's (energy centroids of states with 
given spin / and isospin T) and Et^s (energy centroids of 
states with given isospin T) for proton and neutron sys- 
tems, where one will see that Ei^t^s are approximately 
linear in terms of /(/+ 1) and that Et is precisely linear 
in terms of T{T + 1). In Sec. HI, we present results of 
Ei with random three-body interactions for sd bosons, 
where one sees that regularities of energy centroids of spin 
/ states under random three-body interactions are very 
similar to those under two-body interactions. In Sec. IV, 
we discuss our results of energy centroids with fixed ir- 
reducible representations, where one will see that energy 
centroids with lowest and highest irreducible representa- 
tions carry most of the probability to be the lowest. In 
Sec. V we discuss and summarize the results obtained in 
the present work. In Appendix A we present a few for- 
mulas which are useful in deriving propagation equations 
of this paper. 

In this paper we take the Two-body Random Ensem- 
ble (TBRE) for two-body matrix elements, with the same 




FIG. 1: Probability that an energy centroid Ei^t is the low- 
est. One sees that this probability is large only when I and T 
are close to minimum or maximum. Note that the /max value 
for large T becomes smaller as T increases, and that for 12 
nucleons in the ad shell only / = is possible when T = 6. 



definition given in Ref. j^. For three-body random in- 
teractions for sd boson Hamiltonian, we take the same 
definition given in Ref. p^ : in Appendix B we present 
the definition of three-body Hamiltonian for sd bosons, 
for the sake of convenience. 



II. 



ENERGY CENTROIDS OF SPIN / AND 
ISOSPIN T STATES 



In this Section we investigate regularities of Ej t- Our 



Ej^T values are obtained by using Ei of all systems with 
Np valence protons and Nn valence neutrons under the 
requirement Np + Nn — N . We first obtain number 
of states with fixed / and T, denoted by -D/t, which 
is given by Di{Np = ^i^,7V„ = ^2^) - Di{Np = 
n+2T+2 ^ jy^^ ^ ii^^^^) where Di is the number of 
spin / states with Tz = {Np — Nn)/2. We denote 

n-2T^ 



n-2T\ 



) by eJ^-'^, and Di{Np 



n+2T AT _ 
2 ; * n - 

jjTz-T ^ we can write Ej^t explicitly as follows 



2 )hy Dj''-^ . Using these eJ^-'^ and 



E 



I.T 



Ej^=^ X Dj^-^ 



/D 



IT- 



We have obtained Ej t based on iS/'s of systems with 
TV = 12 (TVp = 6, 7, • • • 12, Np + Nn = 12) and those with 
N ^ 8 {Np = 4, 5, • • • 8, TVp + 7V„ = 8) in the sd shell. 
In Fig. 1 we present out results of V{I,T) which is the 
probability that the lowest energy centroid has spin / and 
isospin T. One easily notices that 'P{I,T) is sizable only 
when both / and T are close to their minimum or maxi- 
mum values. To see this more clearly, we list in Table I 
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TABLE I: J2i T^i^' T) for eight and twelve particles in the sd 
shell. One sees that T) is very small when the value 

of r is not close to Tmin or Tmax. We define ^Y.ii'^^ + 
1)E7^Dit/ {^j{2I + 1)Dit) and V{T) be the probability 
for Et to be the lowest in energy, and note here that ViT) 
does not equal Vjl , T). 



T 
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6 




47.1 
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5.5 


= 8 
2.2 


34.2 








57.9 


0.8 


TV 

0.5 


= 12 
4.9 


0.5 


13.6 


21.8 
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FIG. 2: (-E/T)min versus /(J + 1) values. Here {i5/T)min are 
obtained by averaging the Eit over the cases in which the 
lowest Eit has J = 0, 1 and T = 0. One sees that (i5/T)min 
is proportional to /(/ + 1) approximately, and that (i?/T)min 
with different T values are different: those with larger T being 
systematically higher. 



the values of '^j'P{I^T), where we see that 
is very small when the value of T is not close to its mini- 
mum or maximum value. We note here that ^jV^ItT) 
does not equal V{T), the probability that Et (the energy 
centroid of all states with isospin T) is the lowest. 

The feature that J^i ^(^i T) is large only when T and 
/ are close to their minimum or maximum values is very 
similar to that of 7^(7) discussed in Refs. and more 

generally, to that of 'P(A) discussed in Ref. [13 where A 
denotes irreducible representation of states in interacting 
boson models suggested in Refs. 

As in earlier works, we investigate here the relation 
{Ei)min oc /(/ + 1), where (-E/)inin is obtained by av- 
eraging Ej over the cases of the ensemble in which 
Ej (/ ~ /mill) is the lowest in energy. From Fig. 2, one 
sees that this relation is very robust with inclusion of the 
isospin degree of freedom T, i.e., (i?/,T) min oc /(/+ 1). A 



new and interesting observation here is that the (i?/,T)min 
results, which are obtained by averaging over the case in 
which Eit {I /min, T ~ Tmin) is the lowest energy, can 
be classified according to their T values: the values of 
(£^/,T)min with larger T are systematically higher. 

As we will show later, closer inspection of our 
calculated results confirms that Et = ^iC^I + 
l)Ei^DiT/ (E/(2/ + 1)£'7t) = Eo + CT{T+l) for each 
individual run, which was shown by French many years 
ago Ullll. Dt = X;/(2/+l)£'7T- Note that the {21+1) 
factor is essential for proper definition of fixed-T cen- 
troids. 

Below we discuss fixed-T centroids by propagation 
equations. With nucleons occupying say (ji, ^2, ■ • •) or- 
bits, the spectrum generating algebra is U{N), N — 
2(2ji -f 1) with the factor 2 appearing due to isospin. 
For n nucleons with isospin T , the T quantum number 
labels the irreducible representations (irreps) of SU (2) al- 
gebra that appears in the direct product (space-isospin) 
subalgebra U{N/2) ® SU{2) of UiM). Then the irreps of 
U{Af/2) are completely specified by {n,T). A one plus 
two-body hamiltonian H = h{l) + V(2), which preserves 
angular momentum and isospin is defined by the single 
particle energies (spe) and by the two-body matrix ele- 
ments V/jli = {{kl)Jt I V{2) I iij)Jt), where \{ij)Jt) are 
anti-symmetrized two particle states. With H — h{l) + 
V{2), Et are polynomials in the scalars particle number 
n and T{T + 1]: ^ = ao + ain + a^n^ + az T{T -t- 1); 
see Refs. 0,^^. Solving for the a^'s using Ei^n,T) for 
n < 2, one obtains the following propagation formula 



ET=n{h{\)Y 
' n(n + 2) 
8 

3n(n- 2) 



1) 

2 

r(T + ] 



(l^(2))2." 

(y(2))2^i 



(3) 



with 



From Eq. 
(F(2))2'°} {T, 



i 

* i>j;J 

= H(^(2)>' 



ax 



Ej 



1) 



- T(T+ 1)}. With the inter- 
action matrix elements Kwiw chosen to be zero centered 

*j '■J 

independent Gaussian random variables, ground states 
will have T = or Tmax, with 50% probability for each 
of them. 



III. ENERGY CENTROIDS OF SPIN I STATES 
UNDER RANDOM THREE-BODY 
INTERACTIONS 



The outcome of random three-body interactions for sd 
bosons was first studied by Bijker and Frank in Ref. , 
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FIG. 3; (a) with pure random three-body interactions 

for sd boson systems. One sees that V{I) is large when I 
equals /max and Imin. Vil) also exhibits an odd-even stag- 
gering: ■P(/) is relatively large for odd / values and small 
for even I values, (b) Ei versus /(/ -|- 1). One sees that 
_B/ ~ _Bo + CI{I + 1) for sd bosons with three body random 
interactions. In this figure the total boson number n is 7. 
Similar results are obtained for n = 6 — 20. 



where it was found that the inclusion of random three- 
body interactions does not drastically change the pattern 
of spin I distribution in the ground states, in comparison 
with the results calculated by using random two-body in- 
teractions, if boson number n is much larger than three. 
In this section wc study whether or not the pattern of 
Ei becomes different if one includes random three-body 
interactions. To highlight the feature of Ei with three- 
body part, we use a Hamiltonian with only three-body 
interactions defined in Appendix B. The three-body in- 
teraction parameters are chosen to be random and follow 
the Gaussian distribution. 

Figure 3(a) is a typical example of probability (de- 
noted by V{I)) for Ej to be the lowest with pure ran- 
dom three-body interactions of seven sd boson systems. 
Interestingly but as expected, one sees that V{I) is large 
when / or / ~ /max- One also sees a very appar- 
ent odd-even staggering of V{I) values, i.e., V{I) is large 
when / is odd and relatively smaller when / is even. This 
behavior was also noticed and discussed in Ref. Q when 
only the TBRE was used. 

Figure 3(b) presents (/i^/)min versus /(/ -I- 1) for seven 
sd bosons. A linear correlation between (/?/)inin and 
/(/ -|- 1) can be easily seen. 

A difference between (/?/)min under the TBRE and 
that obtained by using random three-body interactions is 
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FIG. 4: The value of C in the relation {Ei)min ~ C/(/ + 1) 
for sd-boson systems versus boson number n. We see that the 
value of C is approximately proportional to particle number 



found when one studies the correlation (/?/)min — CI{I+ 
1) for sc?-boson systems with different particle number n. 
For the case by using the TBRE, the value of coefficient 
C is not sensitive to particle number n but sensitive to 
the single-particle levels in which the valence particles are 
active; for systems in which one takes random three-body 
interactions, the situation becomes different. In Fig. 4 
we see that the value of C increases with boson number 
n linearly with small fluctuations. 

Now let us investigate the energy centroids of spin / 
states by using propagation equations for random three- 
body interactions. In order to understand the particle 
number dependence seen in Fig. 4, we consider spin / 
centroids Ej generated by random three-body hamilto- 
nians for identical nucleons in a single-j shell. Firstly 
/?/'s correspond to averages over the space defined by 
the irreps n and / of U{2j + 1) and S0{3) respectively 
in U{2j + 1) D 5*0(3). We start with the approximate 
formula (see Eq. (7) in Ref. [13), 



Er = 



(4^g(3))" 

m 



2\n] 



(4) 



which should be valid for j >> n >> 3. In Eq. 10}, 
(//(3))" is the average of H over n particle space and H is 
H with the average part removed (this is made more clear 
ahead). Denoting three particle antisymmetric states by 
\{j)^;al) with a being the extra label required to com- 
pletely specify the states, diagonal 3-particle matrix ele- 
ments of £r(3) are G„,/ = {{if\aJ \ H{3) \ {jf\aJ). 
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To proceed further it is necessary to consider the ten- 
sorial decomposition of and H operators with re- 
sect to U{2j + 1) and the tensors are denoted by v 
[itL ITsI . a general fc-body operator wiU have v = Q, 
1, . . ., k parts. However for a single j shell the v ~ 1 
part will be zero. Thus P = (/2)>.=o _^ ^2y=2 ^nd 
H{i) = iJ'^=°(3) + H''=^{3) + H''=^{3). The and 
(I^)" — are generated by the v = Q parts, and 



7i(2j + l-n)(j + l)~ + 



(5) 



Note that H = H - H''=°, and (/^^(S))" = 
{{PY='^H''=^{3))'^ = ((/2)'^=2^(3))". It should be rec- 
ognized [13 that 77-^=2(3) = _ 2)^^=2 (2), where F is 
a two-body operator with rank v = 2 and n is number 
operator. For our purpose, it is not necessary to know 
the exact form of the F operator. The propagation equa- 
tion for the n particle average {{PY^'^ F'^^'^ {2))'^ is well 
known (see for example Eq. (A2) of Ref. Using 
this we obtain 



n{n - l){n - 2){2j + 1 - n){2j - n) 



2(2j-l)(2j-2) 



(6) 



Putting n = 3 on both sides of Eq. ®, the average 
involving F can be eliminated and this gives 



(2j + 1 ■ 



_ n(n-l)(n-2) 
" 6(2j-2)(2j-3) 
n){2j-n){{PY^^H{-i)f 



(7) 



As given in Ref. T^, (I^) 



u=2 



K2j-fl-n)(j + l). 



Combining Eq. with Eqs. Q and give the final 
formula for E(^n j^ with the /(/ + 1) term carrying linear 
n dependence 



E^ 



{n.I) 

3n 



Ea AAJ + 1) - 3j-(j + l)}G„j(2J + 1) 



[j(j + l(2j + l)]2(2j + l) 
X /(/ + !) , 



(8) 



where Eq is a constant determined by Gaj, n and j. 
The form of Eq is given in Appendix A. The second term 
in Eq. (jS)) clearly gives a linear n-dependence for C. 
This dependence, as seen from the unitary decomposi- 
tion, comes from the v — 2 part of H{3) which in turn is 
responsible for the /(/ -I- 1) term. 

Similarly, one investigates E(^n,i) versus /(/ + 1) for 
bosons with spin I. Eq. Q remains the same but 



(/2)" ^ -n{2l + l + n)l 



n{n - l)(n - 2)(2/ + 1 + n){2l 2 + n) 



{21 + 1){21 + 2)(2l + 3)(2l + A){21 + 5) 
+ 1) - 3l{l + 2)}G„l(2L + 1) . 



(9) 



Then one obtains 



E, 



{n.I) 



E'q 



P{2l + l){2l + 2)---i2l + 5) 
J2 [L{L + 1) - 3l{l + 2)] G«i X /(/ + 1) , (10) 

where GaL are three-body matrix elements for bosons 
with spin I. For d bosons there are five three-body ma- 
trices with L = 0,2,3,4 and 6, respectively. By using 
Eq. (|10() , one obtains the value of coefficient C in the re- 
lation (£'/)min — C/(J -|- 1) for d boson systems: 



6n 



TrP{2l + l)(2l + 2)---(2l + 5) 
((L(L + 1) - 3/(/ + 2))(2i + 1))' 



1680 



n ~ 0.03331n 



(11) 



where the constant \/ - comes from the integral 



J xexp{ — -) dx 



The value of C obtained by our numerical experiments 
(1000 runs of the ensemble with random three-body in- 
teractions for d boson systems of n = 6-20) is 0.03499n. 
The value of C is therefore reasonably reproduced by the 
propagation equation (which predicts C = 0.03331ri, as 
discussed above). 



IV. ENERGY CENTROIDS WITH FIXED 
IRREDUCIBLE REPRESENTATION 

In this section we consider two interesting examples: 
(i) energy centroids with fixed irreps of the SUsd{i) ® 
SUr ,f(3) limit of sdpflBM which was mentioned in Ref. 
|10| : (ii) energy centroids with fixed Wigner's spin-isospin 
supermultiplet SU{4) irreps for {2sld) shell nuclei. 

Let us first consider energy centroids with fixed ir- 
reps [nsdiKd^J'sd) ■ npf{\pfiipf)] of [Usd[&) 3 SUsdii)] © 
[Upf{lQ) D SUpf{3)\ subalgebra of the spectrum gen- 
crating algebra (SGA) [/sdp/ (16) of sdp/IBM; see Refs. 
[Ii for details of the SU{3) hmit of sdp/IBM. Given a 
one plus two-body sdpf hamiltonian, with boson number 
n = risd + npf where Usd and Upf are number of bosons 
in sd and pf orbits, the propagation equation for energy 
centroids can be written as 



E.. 
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= ao + aiHsd + a2npf + 03 ^ "^^'^ j + 04 ^ ^'^■'^ j 
+a5nsdnpf + 066*2 (Asd^sd) + a7C'2(Ap/^p/) , (12) 

where C2({/}) in this section denotes the eigenvalue of 
the quadratic Casimir invariant of a given irrep {/}. For 
irreps (A^) of SU(3) it is given by C2(A^) = X"^ + fi'^ + 
Xfi + S{X + fi). The final propagation equation can be 
obtained by solving for the 's in terms of the centroids 
with n < 2. We give this equation in the Appendix A. 

In order to calculate energy centroids using Eq. ((T^ . 
the reductions Usd — ^ (AsdMsd) and Upf (Ap//ip/) are 
needed. For sd bosons the reductions are well known 
[1^ ITsI and for the pf system the reductions are ob- 
tained using the method given in Ref. 20]. Following 
the results for sdlBM-T and sdlBM-ST energy centroids 
given in Ref. 10], we consider the basic energy cen- 
troids with n < 2 as indepen- 
dent zero centered (with unit variance) Gaussian random 
variables, instead of considering the single-particle en- 
ergies and two-body matrix elements in sdpf space as 
random variables. 

With 1000 samples, the probabilities for the energy 
centroid with a given [usdiXsdfJ'sd) ■ n-pf{Xpffipf)] irrep 
to be lowest are calculated for boson numbers n = 8, 9 
and 10 and the results are shown in Fig. 5. Firstly it 
is seen that the irreps with the lowest and highest Usd 
carry most of the probability, about 84%. For each of 
the other Ugd the probability is 1 — 3%. Moreover for 
the lowest and highest Usd, the probability splits into 
the lowest and highest SU{3) irreps. For Ugd = ob- 
viously {Xsd^J'sd) — (00) and the probability for highest 
(according to the eigenvalue of the SUpf{3) quadratic 
Casimir invariant) SUpf{3) irreps (3n,0) is ^ 24% and 
for the lowest irreps it is ^ 19%. The lowest irreps for 
Upf = n = 8, 9 and 10 are [Xpffipj) = (00), (30) + (03) 
and (00), respectively. Similarly, for Usd = n, the high- 
est SUsd{-i) irreps are (2n, 0) with probability 24% and 
the lowest irreps are (02), (00) and (20) respectively for 
nsd — n — 8, 9 and 10 with probability ~ 17%. Thus 
the results in Fig. 5 show that the energy centroids with 
the lowest and highest [nsd{Xsd^J'sd) ■ "'p/(Ap//ip/)] irreps 
carry most of the probability just as with / and IT en- 
ergy centroids considered in this paper and in many other 
examples considered in Refs. [7-12]. 

Now let us come to the second example, energy cen- 
troids with fixed SU{4) - (ST) irreps for the (2sld) shell 
model. For {2sld) shell nuclei, spin-ispopin supermulti- 
plet SU{4) algebra appears in the direct product subal- 
gebra U{6) ^ SU{4) of [/(24) SGA. 

We first note that U{6) generates the orbital part and 
SU{4:) generates spin-isospin (ST) quantum numbers via 
SUsri'i) 3 SUs{2) (g) SUt{2). For a given number of nu- 
cleons n, the allowed U{4) irreps are {/} = {/i, /2, /a, fi\ 
with h>f2>h>h> 0, /i < 6 and f^ + f^ + h+U^^ 
n and the C/(6) irreps, by direct product nature, are {/}, 
the transpose of {/}. It is important to note that the 
equivalent SU{A) irreps are {/i - /4,/2 - h, h - U}- 
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[n(2n,0):0(00)] — ■ y 
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FIG. 5: Probabilities for sdpf IBM centroid energies with 
fixed-[nsdiXsdfJ'sd) ■ npf{\pfiipf)] irreps to be lowest in en- 
ergy vs for systems of n = 8, 9 and 10 bosons. For 
riad = (ripf = n) and Usd ~ n (upf = 0) the probabilities 
for (XpfUpf) [{XsdfJ-sd) = (00)] and (AsdMsd) [(-^p/Mp/) = (00)] 
are shown in the figure. However for < rigd < n, the 
probabilities shown are the sum of the probabilities for the 
[nsd{X3d^J-sd) ■ Wp/(Ap//ip/)] irreps; for each Usd and Upf = 
n — Usd there are 3 — 4 SU{3) irreps but they are not shown 
in the figure to avoid clustering of too many points. Filled 
circles, squares and stars are joined by lines to guide the eyes. 
See text for further details. 



With these, from now on we will use /7(4) and the ir- 
reps {/}. It is well known that a totally symmetric U{A) 
irrep {A} ^ {ST) = (|, A), (A - 1, A _ 1), . . (QO) 

or (^, |). Using this result and expanding a given ?7(4) 
irrep into totally symmetric U (4) irreps will give eas- 
ily {/} i^^r) reductions. Just as the fixed-T energy 
centroids propagate, the fixed {f}{ST) energy centroids 
^{f}{ST) for a one plus two-body hamiltonian propagate 
as the available scalars of maximum body rank 2 are 1, 
n, n^, C2(C/(4)), 5*2 and T^ and the centroids for n <2 
are also six in number. The propagation equation, with 
C2{{f}) ^E^f! + 3/1 + /2 - A- 3/4 for «7(4) irrep 
{/}, was first discussed in Ref. |2j] and we present it in 
Appendix A, for the sake of convenience. 

Just as in the sdpf example, we consider the basic 
energy centroids (H)^^^^^'^^ with n < 2 as independent 
zero centered (with unit variance) Gaussian random vari- 
ables, instead of using ci and V^^li as random variables, 
and study the {f}{ST) structure of the ground states. 

Using 1000 samples, the probability for a given fixed- 
{f}{ST) energy centroid to be lowest in energy is calcu- 
lated and the results are shown in Fig. 6 for n = 8, 9, 
10 and 12. The probabilities split into three ?7(4) irreps 
(other irreps carry < 1% probability) for n = 8, 9 and 
10, and the corresponding (ST) values are as shown in 
Fig. 6. Energy centroids with the lowest and highest 
U (4) irreps carry ~ 25% and ~ 40%, respectively. The 
lowest irreps are {2^}, {32'^} and {3^2^} respectively for 
n = 8, 9 and 10 and the highest irreps are {6,n — 6}. 
The third irreps {4^}, {54} and {5^}, with probability 
~ 32% for 71 = 8, 9 and 10, are those that carry S — n/2 
or T = n/2. Besides these, for n ^ 10 the irrep [33l](00) 
carries 3.7% probability. For the mid-shell example with 
n = 12 the probabilities split into the lowest {3^} and 
highest {6^} irreps with ~ 25% and ^ 75%, respectively. 
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• n=8 : 
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FIG. 6: Probabilities for (2sld) shell centroid energies with 
fixed {f}{ST) irreps to be lowest in energy vs C2({/}). Re- 
sults are shown for nucleon numbers n = 8, 9, 10 and 12. The 
irreps carrying probability less than 1% are not shown in the 
figure. The 17(4) irreps {/} for the results in the figure are 
given in the text. The corresponding (ST) values are shown 
in the figure as {2S2T). Filled circles, filled squares, open 
circles and open squares are joined by lines to guide the eyes. 



The lowest irrep supports only (ST) = (00) and the prob- 
ability for the highest irrep spUts into ~ 13% and ~ 62% 
for (ST) = (00) and (12,0) -I- (0,12). A very important 
observation from Fig. 6 is that the probability for the en- 
ergy centroid with lowest C/(4) irrep to be lowest is only 
^ 25% and it should be noted that the corresponding 
SJ7(4) irreps are {0}(00), {l}^) and {1^}(10) + (01) 
for n — 4fc, 4fc + 1 and 4fc4-2, respectively, with k being a 
positive integer. In fact as discussed in Ref. |2^, realistic 
interactions give ground state wavefunctions having over- 
lap of ^ 90% with these irreps, i.e. very high probability 
for a-cluster structure. However detailed calculations in 
Ref. I23 showed that random interactions give a very 
small probability for a clustering. The same result has 
been brought out in a simple and easy manner by the 
{/}(ST) energy centroids. 

Finally we point out that the present study extends 
easily to (2pl/)-shcll nuclei by changing the restriction 
/i < 6 to /i < 10 in enumerating the U{4) irreps. 



dom three-body interactions. Results of T'(/), the prob- 
ability for Ej to be the lowest in energy, and (-E/)inin 
(and (£^/)max), obtained by using the TBRE, are also 
applicable to those by using random three-body interac- 
tions, except that the value of coefficient C in the relation 
(i?/)niin — CI{I + 1) increases linearly with boson num- 
ber n. 

Third, such regular patterns can be generalized to en- 
ergy centroids with given irreducible representations of 
groups for boson systems as well as fermion systems. 
We consider energy centroids with fixed irreps of the 
SUsdi^) ® SUpf{3) limit of sdpflBM, and energy cen- 
troids with fixed Wigner's spin-isospin supermultiplet 
SU{A) irreps for {2sld) shell nuclei. We see that the 
lowest and highest irreps {/} carry most of the cases 
that Ej is the lowest in energy, and the energy centroids 
propagate via quadratic Casimir invariants. 

The above results, such as energy centroids with fixed 
T value, energy centroids of spin / states under random 
three-body interactions, energy centroids of fixed irreps 
of the IBM models and shell models, are discussed by 
using propagation equations. 

Our results suggest that behavior of energy centroids 
discussed in Refs. [7-12] is a very robust feature for quan- 
tum many-body systems interacting by random forces. 
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V. DISCUSSION AND SUMMARY 



In this paper we studied the behavior of energy cen- 
troids in the presence of random interactions. First we 
show that results of energy centroids are robust regard- 
less of inclusion of isospin, based on numerical calcula- 
tions. We find that (i?/T)inin (and (i?/T)max) values can 
be classified according to their T values. The simple re- 
lation, Et = Eq + CT{T + 1) for each individual run, is 
confirmed and discussed. We see Ejt with both neutrons 
and protons shows a similar pattern to that of Ej with 
only identical particles discussed in earlier works. 

Second, we find in this paper that for sd boson sys- 
tems the feature of i?/'s is robust with inclusion of ran- 



8 



[1] C. W. Johnson, G. F. Bertsch, and D. J. Dean, Phys. 

Rev. Lett. 80, 2749(1998). 
[2] Y. M. Zhao, A. Arima, and N. Yoshinaga, Phys. Rep. 

400, 1 (2004). 

[3] V. Zelevinsky and A. Volya, Phys. Rep. 391, 311 (2004). 
[4] T. Papenbrock and H. A. Weidenmueller, Phys. Rev. 

Lett. 93, 132503 (2004); AIP Conf. Proc. 777, 140 

(2005). 

[5] N. Yoshinaga, A. Arima, and Y. M. Zhao, "Lowest 
bound energies for random interactions and the origin 
of spin zero ground state dominance", preprint (to be 
published). 

[6] T. Otsuka and N. Shimizu, AIP Conf. Proc. 726, 43 
(2004). 

[7] A. Arima, N. Yoshinaga, and Y. M. Zhao, Eur. Phys. J. 

A 13, 105 (2002); N. Yoshinaga, A. Arima, and Y. M. 

Zhao, J. Phys. A 35, 8575 (2002). 
[8] Y. M. Zhao, A. Arima, and N. Yoshinaga, Phys. Rev. C 

66, 064323(2002). 
[9] Y. M. Zhao, A. Arima, and K. Ogawa, Phys. Rev. C 71, 

017304 (2005). 
[10] V. K. B. Kota, Phys. Rev. C 71, 041304 (2005). 
[11] D. MulhaU, A. Volya, and V. Zelevinsky, Phys. Rev. Lett. 

85, 4016(2000). 

[12] V. K. B. Kota and K. Kar, Phys. Rev. E 65, 026130 
(2002). 

[13] V. Velazquez and A. P. Zuker, Phys. Rev. Lett. 88, 
072502 (2002). 

[14] A. Arima and F. lachello, Ann. Phys. 99, 253 (1976); 
ibid. Ill, 209(1978); ibid. 123, 468 (1979). 

[15] F. lachello and A. Arima, The Interacting Boson Model 
(Cambridge University Press, England, 1987). 

[16] R. Bijker and A. Frank, Phys. Rev. C 62, 014303 (2000). 

[17] H. Banerjee and J.B. French, Phys. Lett. 23, 245 
(1966); J.B. French, in Isospin in Nuclear Physics, edited 
by D.H.Wilkinson (North Holland, Amsterdam, 1969), 
p.259; F.S. Chang, J.B. French, and T.H. Thio, Ann. 
Phys. (N.Y.) 66, 137 (1971). 

[18] J.C. Parikh, Group Symmetries in Nuclear Structure 
(Plenum, New York, 1978); S.S.M. Wong, Nuclear Statis- 
tical Spectroscopy (Oxford University Press, New York, 
1986). 

[19] J. Engel and F. lachello, Phys. Rev. Lett. 54, 1126 
(1985); H.Y. Ji, G.L. Long, E.G. Zhao and S.W. Xu, 
Nucl. Phys. A 658, 197 (1999). 

[20] V.K.B. Kota, J. of Phys. A 10, L39 (1977); V.K.B. Kota, 
Physical Research Laboratory Technical Report PRL- 
TN-97-78 (1978); Mathematics of Computation 39, 302 
(1982). 

[21] R. U. Haq and J. C. Parikh, Nucl. Phys. A220, 349 
(1974). 

[22] Y.M. Zhao, A. Arima, N. Shimizu, K. Ogawa, N. Yoshi- 
naga, and O. Scholten, Phys. Rev. C 70, 054322 (2004). 



Appendix A Useful formulas in deriving propagation 
equations 

First we present the detailed result of Eq in Eq. (9) of 



Sec. in. 



En = 



(2j 



a, J 

3n2 



2j(j + l)(2j + l)3 
X ^{ J( J + 1) - 3i(j + l)}G„j(2 J + 1) 



a, J 



Second, we present the propagation equation of energy 
centroids with fixed irreps [nsdi^sd^J'sd) ■ ''T'pfi^pff^pf)] of 
[Usd{6) D SUsdm © [C/p/(10) D SUpf{3)] subalgebra 
of the spectrum generating algebra (SGA) C/sdp/(16) of 
sdpf IBM. This is obtained by solving a^'s in Ea. (|12(l by 
centroids with n < 2. 



E, 



= [1 - Usd - ripf 



nsd 
2 



2 



nsdnpf 



Eo 



(00):0(00) 



nsd-2 { ^^'^ I - Usdripf 



El 



(20):0(00) 



sd'^pf 



E, 



0(00):1(30) 



-71s(jnp/i?l(20):l(30) 

5 / risd 



g""'" ' 9 I 2 

r'^+g[ 2 
3 2 / Upf 
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C2{^sdfJ-sd) 



E- 



2(40):0(00) 



E. 



2(02):0(00) 



Upf 



+ ■^C'^i^pff^pf) 
■^C2{XpfHpf) 



E, 



0(00):2(60) 



E, 



0(00):2(22) ■ 

(13) 



Last, we give the propagation equation for energy cen- 
troids with fixed spin-isospin SU(4) irreps for the sd shell 
nuclei. This was first discussed in Ref. 12111. 



E. 



{f}{ST) 



(l - ^" + l^^^ ^{o}(oo) + (2n - ^{i}(ii) 

^ln+^n^ + ^C,{{f}) + lsiS + l) 
1 



{2}(11) 



-T(r + i) 



-ln+lc2i{f})-\siS + l) 
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^-r(T + i) 



E 



{2}(00) 



+ 



ln+ln'^^C2{{f}) + \s{S+l) 



-T{T+1) 
4 



ln+ln^-lc2i{f})-\siS + l) 



1 



T(r + 1) 



(14) 



Appendix B Definition of three-body interactions of 
sd bosons 

In this Appendix we present the definition of the three- 
body Haniihonian of sd boson systems discussed in Sec. 
III. We note that the same definition was taken in Ref. 
^16] by Bijker and Frank. Discussions of three-body in- 
teractions in the interacting boson model can be found 
in Ref. [13. 



Our three-body Hamiltonian of sd bosons are given by 



^3= E E^^. 

L=0,2,3,4,6 i<j 



PL ■ Pu 



where 





1 

^ 6 


{sKsls^f^ 




1 

^ 2 


(stdtrft)(o) 




1 

^ 6 


{duu^f^ 




1 

^ 6 


(ststdt)(2) 


p^ 


1 

~ 2 






1 

~ 6 


(dtdtdt)(2) 


Pi 


1 

^ 6 


(dtdtdt)(3) 


pt 


1 

^ 2 




p^ 


1 

~ 6 


{dUU^f^ 


pt 


1 

^ 6 


(dtdtrft)(6) 



,(15) 



The coefficients ^l^j are random and follow the Gaussian 
distribution with their widths given by 



(a.,?L',.,> -(1 + '5ll'<^,.'%')/2 



(17) 



